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Abstract 



We show that dark matter could be made of massive gauge bosons whose stability 
doesn't require to impose by hand any discrete or global symmetry. Stability of gauge 
bosons can be guaranteed by the custodial symmetry associated to the gauge symmetry 
V.O an d particle content of the model. The particle content we consider to this end is based 

on a hidden sector made of a vector multiplet associated to a non-abelian gauge group 
and of a scalar multiplet charged under this gauge group. The hidden sector interacts 
with the Standard Model particles through the Higgs portal quartic scalar interaction in 
Qh such a way that the gauge bosons behave as thermal WIMPS. This can lead easily to the 

observed dark matter relic density in agreement with the other various constraints, and 
can be tested experimentally in a large fraction of the parameter space. In this model 
the dark matter direct detection rate and the annihilation cross section can decouple if 
the Higgs portal interaction is weak. 

> 



1 Introduction 



In the Standard Model (SM) the best known sector is the gauge sector. Besides Lorentz 
symmetry and CPT symmetry, all fundamental symmetries of the SM are gauge symmetries. 
There are also accidental symmetries, in particular B — L number conservation (or B and L 
separately neglecting highly suppressed instanton effects), but also this one is closely related 
to the gauge symmetries. It results from the particle content of the SM and the gauge sym- 
metry charges assigned to it. There are no ad- hoc discrete or global symmetries. This leads 
to a certain number of stable particles, the photon (stable due a conserved gauge symme- 
try which makes it massless), the electron (stable because it is the lightest particle charged 
under this conserved gauge symmetry), the lightest neutrino (due to Lorentz invariance and 
the fact that it is the lightest fermion) and the proton (from the accidental baryon number 
conservation). Any new physics model must not spoil the stability of these particles with 
huge accuracy. 

In the following we assume that this gauge symmetry reason for having stable particles 
also holds for the dark matter (DM) particle. Of course there is no mandatory reason why 
this must be the case for DM too, but the fact that this holds for the SM is puzzling enough 
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to investigate this possibility. Along this line we consider a particularly simple model where 
DM is made of a multiplet of vector particles communicating with the SM through the Higgs 
portal, and show it is a perfectly viable DM candidate^ In the process we present a new 
model with spin-1 DM. Most DM models involve particles of spin (scalar singlet [H [2j G], 
axion pQ, inert doublet OE], high SU(2)i scalar multiplets [6],...), spin 1/2 (neutralino [7], 
axino [8j, high SU(2)l fermion multiplets [6], singlet fermion [9l[in|,....) or s Pi n 3/2 (gravitino 
particles. To our knowledge models of spin-1 DM which have been proposed involve 
either more than 4 dimensions [12] or an explicit discrete symmetry (little Higgs DM with 
T-parity p3]). 



2 Model 

If there were no fermions in the SM, and no mixing between SU(2)l and f7(l)y, the W 
and Z bosons would be degenerate in mass and stable due to the custodial symmetry of the 
scalar and SU(2)l gauge SM sectors. As, fortunately, there are fermions (and SU{2)l x U{\) 
mixing), the W and Z are not stable, but this possibility to have stable gauge bosons could 
hold in a hidden sector coupling to the SM through the Brout-Englert-Higgs (Higgs for short) 
portal quartic scalar interaction. 

The scenario is based on two simple assumptions. First we assume the existence of a A'^ 
gauge multiplet associated to a non-abelian gauge symmetry, G' . Under this gauge group all 
SM particles are singlets. Second this non-abelian field couples to the SM only through a 
complex scalar Higgs portal field, eft, which is singlet of the SM but charged under G' . Mixing 
of A'^ with the SM gauge bosons (through F^Fy kinetic mixing with the hypercharge gauge 
field) is forbidden by the non-abelian character of the extra gauge symmetry. For instance 
we consider the simple case G' = SU(2), which we denote SU(2)hs, and we take ^ to be a 
complex doublet of this gauge symmetry. No discrete symmetry are assumed at any level. 
The most general Lagrangian one can write with these simple assumptions is 

£ = C SM _ 1^ . p,^ + (iV) t ( £,/Y) _ Am 0W# " 4<P<t> ~ A^^V) 2 , (1) 

with D^cj) = d^<j) — i^r ■ A'V 1 . In the SM lagrangian we define the Higgs potential notations 
as: C SM 3 -n 2 H^H - X(H^H) 2 with H = (H+, H°). 

If SU(2)hs is spontaneously broken (i.e. < 0), writing as usual, = exp(ir ■ £/i></,) • 
(0, -^[v<f> + v']) 7 \ an d gauge rotating away the £ part to absorb it in A^ = UA'^U^ 1 — 
[3 M C/]{7 _1 with U = exp(—ir ■ C/ V <f>) we get 
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C = C S M-\F^.F^ + ^(g^) 2 A^A^ + ^glA^A^' 2 + ^glv^.A^' 
which gives = g ( j ) v^/2 and m 2 , = —2/^. 



1 Our perspective is purely phenomenological, and therefore different from the top-down one of particular 
grand-unified supersymmetric models which can achieve the same goal, i.e. obtaining R-parity conservation 
as a remnant symmetry of a spontaneously broken gauge group, e.g. gauged U(1)b-l |14j . 
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The Lagrangian above has an important property: since the scalar field is in the fun- 
damental representation of the gauge group it displays a custodial symmetry, 50(3) in the 
^123 component space. As a result the 3 Af components are degenerate in mass and are 
stable. Since all other particles (as rf) are SO(3) singlets, any decay of the gauge bosons is 
forbidden by the custodial symmetry. In practice, this stability can be also seen from the 
facts that: 



• All interactions from the scalar kinetic term involve the A^ field in pairs. The dangerous 
d^cf)A^ terms which would make the gauge boson instable, disappear from the absorption 
of the Goldstone bosons by the gauge bosons, as in the SM. These scalar kinetic term 
interactions do not mix the various Af SU (2) l components. 

• In the pure gauge sector there are quadrilinear and trilinear interactions but they do not 
cause any decay of the gauge boson fields. The quadrilinear terms because they involve 
A^ components in pairs. The trilinear terms because they involve three different Af 
components - they are of the form eijkd fl Ai u (Aj A v k — A v - A^) - while all other interactions, 
which would allow the vector bosons to decay (to scalars), involve two same components. 

This custodial symmetry structure is closely associated to the particle content of the 
model. For instance, light fermions charged under G' or higher scalar multiplets of G' would 
make the gauge bosons unstable, see section 8. But in a similar way could we think about 
new particles destabilizing the proton which nevertheless turns out to be stable with huge 
accuracy. 

In order to get the mass spectrum of this model it is necessary to minimize the potential 
not only along the (j) direction as above but also along the H direction. Writing H = 
exp(ir ■ (/v)(0, ^[v + h']) T , with v = 246 GeV, going to the unitary gauge and imposing 
dV/dh = dV/dr] = we get: 

2 + §<W| 

v = ~Y\ (3) 



^ AA,A — tAJL 



4 "to 



This leads to a non-diagonal mass matrix in the (h' , rf) basis 



ml, m 2 h . 



m 2 = y y , (5) 



with 



m h> v > m v > 



-2^ 2 AA^ - |AA m /4 + §AA m /4 
mh ' = XX, ^ > ( 6 ) 

*A(j> ~ 4 A m 



2 ~~ 2^/?AA0 — ^XfpXm^i 2 + |A^A m ^ 2 
"V = AA, *A 2 ' () 

m h'r,' = ^mWtj,. (8) 
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Figure 1: Annihilation diagrams with no DM particle in the final state 



A, A k A ' 

A > ^^' h Aj 

Figure 2: Annihilation diagrams with one DM particle in the final state 

Diagonalizing the mass matrix with 

( h ' \ _ ( cos P s™/ 3 \ ( h 
\ rf J \ — sin (3 cos /3 J \ rj 

where tan 2/3 = 2m^, -/(mj, — iri^,), one obtains the following Lagrangian in the physical 
state basis (omitting all SM terms which are unaffected by the h'-rf mixing) 

+ A^A" Irfi] 2 + 4h 2 + 4 vV h + 2v^V + 2i>^fc] 

+ (2W+W-» + 1 Z^Z^rf 2 + K h h 2 + K hvV h + 2v^r t + 2v£ h h] 
p cos^ Vw 

- X v r] 4 - X h h A - Xxtfh 2 - X 2 h 3 r] - X 3 hr] 3 - p v r, 3 - p h h 3 - Pl r] 2 h - p 2 h 2 i] . (10) 

The various parameters of this Lagrangian are given in appendix A in terms of the input 
parameters of Eq. ([!]). Note that, beside the SM parameters, this model has only 4 free 
parameters what makes it potentially well testable: the four input parameters g^, X^, p^ and 
X m or equivalently the masses, symmetry breaking scale and gauge coupling, iua, m v , and 
9<t>- 

3 Relic density 

There is a number of tree level annihilation processes which determine the vector multiplet 
relic density through their freezeout. They can be classified in two classes. First the annihi- 
lations with no trilinear gauge interactions: AiAi — > rjrj, rjh, hh via a direct quartic coupling 
or via an intermediate h or rj; and AiAi — > //, W + W~ , Z°Z° via a h or 77 in the s-channel, 
Fig. 1. Second the channels with one trilinear gauge coupling, Fig. 2. These channels have 
the peculiar property (impossible with ordinary models based on a Z 2 symmetry) to involve 
one DM particle in the final state by reducing the number of DM particles from 2 units to one 
unit. In the following, at the exploratory level of this article, we calculate the relic density 





numerically, limiting ourselves to the diagrams of Fig. 1 which in most of the parameter space 
are expected to give a relatively good approximation^] 

In the following we will require that the relic density obtained is within the 3a range 
0.091 < nh 2 < 0.129 [I?]. 



4 Direct detection 

At tree level a vector DM particle can collide elastically a nucleon either through h exchange 
or via 7] exchange, which results in a spin independent cross section 

a SI (NA -> NA) = -^sin 2 2fim% ^~ ^ 4 . (") 
b47r ^ f z m^m^ m A 

with m r = m^iTiA/ '{tun + mA) the reduced mass and mjy the nucleon mass. / parametrizes 
the Higgs nucleon coupling, frriN = (iV| J2q m qQl\N) = 9hnnv. We take the value / = 0.3. 
For rriA » mN numerically one gets: 

v ; M V0.3/ V0.5/ V500GeV/ V ra A / V / 

(12) 

with i? = (m 2 — mfj 2 /m^ which is unity for m v » m^. This has to be compared with 
the present experimental upper bound on this cross section [18] which, for example around 
ttidm = 100 GeV, is of order 10 -44 cm 2 . Therefore for not too small /? mixing angle the 
expected signal can be easily of order the present sensitivity or even exceed it by few orders of 
magnitude. For small mixing angle it can be well below it, independently of the relic density 
constraint (see below). Direct detection constraints are therefore already relevant to exclude 
a part of the parameter space. 



5 Electroweak precision measurement constraints 



The main constraint comes from the contribution of the 77 scalar to the T parameter which 
is the same as the one of a scalar singlet mixing with the Higgs boson [191 [2] : 
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2 We leave to a subsequent publication [15] the analytic determination of all the cross sections as well as 
the incorporation of the diagrams of Fig. 2. These diagrams are of same order in than the first and third 
diagrams of Fig. 1 and do not cause any particular large effects in the Boltzman equations. This means that 
in order to get the observed relic density their contribution can be compensated by a moderate decrease of 
the coupling (and therefore a moderate increase of for tubm fixed). Similarly the incorporation of 
these scatterings will not spoil the two operative regime picture explained below. To incorporate these effects 
necessitates a modification [15] of the source code of the program MicroMega2.0 [16] we have used. 
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where s 2 = sin 2 9y/{Mz) gives the weak mixing angle in the MS scheme at the scale \i = Mz, 
c 2 = Myy/Mg, and rriH is the reference value of the Higgs boson mass in the SM. The 
phenomenology of Eq. (13) has been studied in Ref. |19| (see Figs. 9-10 of this reference 
where similar expressions for the less important S and U parameters can also be found). For 
rrih = rrirj the constraints on rrih and m v are the same as the ones on the Higgs boson mass in 
the SM: rrih < 154 GeV at 95% confidence level |20j. For maximal mixing if one scalar mass is 
below this value the other one can be larger but should remain under ~ 250 GeV. For smaller 
mixing angle the rf dominated scalar can be easily much heavier but the one predominantly 
composed of b! must remain low. In the following we will limit ourselves to require that 
T — Tsm is within the conservative range —0.27—1-0.05 |21j (taking mn = 114.4 GeV in 
Eq. (13)). For < 114.4 GeV we will also require that the rj — >■ // branching ratio (or 
equivalently sin 2 f3) is below the upper bound from direct search at LEP, Fig. 10 of Ref. |22j . 



6 Results 

6.1 Small Higgs portal coupling 

If A m is small, but large enough to thermalize the rj and A^s with the SM thermal bath prior 
to DM freeze out, (for instance for example for rriA ~ 100-1000 GeV, within the range say 
X m ~ 10~ 8 — 10~ 3 ), the h'-rj' mixing angle /5 is small (except for rrih' — rn v > in case the mixing 
can be large independently of the size of the off-diagonal term in Eq. ([5|). Except for this 
case, and except also for rrih ~ 2rriA or rrih ~ 2rriA (where the h or r\ exchange diagrams can 
be resonantly enhanced and therefore be relevant even for small f3 angle), the only relevant 
process for the relic density in Fig. 1 is the A{Ai — > rjr] process. In this case this process 
depends only on g^, and A^. If A<^ is in addition small, only the first and third diagrams of 
Fig. 1 remains and the annihilation cross section depends only on g^ and rriA (or equivalently 
v^). The dependence in the small rj mass, ~ ■^2\ ( f ) v ( f > can be neglected. This leads to 
a cross section proportional to g^/m A , that is say proportional to g 2 /^ 2 , and therefore to a 
linear correlation between and given in Fig. 3. a (red dots for instance). Assuming a 
perturbative g^ coupling leads therefore to an upper bound which is of order few tens of TeV 
depending on the perturbativity condition considered, we get for example rriA < 25 TeV if 
g^ < 4-7T (in agreement with the unitarity bound which holds for any thermal particle whose 
relic density results from the freezeout of its annihilation |23j). The corresponding values of 
rriri versus rriA are given in Fig. 3.b. scales like s/rriA due to the linear correlation between 
g<f, and v$. Both rriA and m r; can be as small as 1 GeV or even much lessj^] For low h-rj 
mixing such low values of the masses are allowed by electroweak data as well as by direct 
search limits from LEP. The constraints on the Higgs mass, to a good approximation, are 
the same as in the SM. As for the direct detection rate, it is proportional to the small Higgs 
portal interaction X m , and consequently essentially decouples from the relic density which, as 
shown above, is essentially determined by the pure hidden sector parameters independently 

3 We don't plot what happens below 1 GeV because the code Micromega we used doesn't allow us to go 
much below this scale. But the linear correlation above between and v$ clearly holds to much lower values 
as long as is small enough. The masses must be nevertheless above the ~ MeV scale due to BBN constraints 
|24| . Whether this could lead to an explanation of the 511 KeV gamma rays from the galactic center observed 
by Integral [25] would be worth to be analyzed (through DMDM — > rjij annihilation followed by rj decays to 
e+e"). 
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Small A,„ regime Small A,„ regime 




g* m A (GeV) 



Figure 3: vs and uia vs m v leading to 0.091 < Qh 2 < 0.129, for 10~' < A m < 10~ 3 , 
rrih = 120 GeV and various values of A^: = 10 -4 (red), A^ = 10~ 3 (orange), A^ = 10~ 2 
(green) and A^ = 10 _1 (blue) (from left to right and top to bottom respectively). One also 
recognizes the uia = g^v^jl ~ m/j/2 resonant case curve. 



of A m As a result, for 10 7 < A m < 10 3 , it is e.g. few orders of magnitudes below the 
present upper bounds on the direct detection rates [18]. This decoupling allows to avoid 
the tension which exists between both constraints in various models. Note however that for 
rrirj « rrih (i.e. small A^), the cross section of Eq. 

so that even for small X m we can get a large direct detection cross section if A^ is even more 
suppressed. For example for A^ as small as 10~ 4 and A m = 10~ 3 we get a cross section which 
can be as high as 10 -43 cm 2 independently of tua- 

For larger value of A^ the second diagram of Fig. 1 also becomes important which modifies 
the correlation between and v^, Fig. 3. Larger values of g§ are necessary for small v^. 
For large value of A^ in order to have enough suppression of the annihilation cross section 
one needs either m v close to tua or uia in the multi-TeV range (to benefit from the l/m A 
asymptotic behavior of the cross sections), or one must have one of the 2 resonances above 
effective. 



11 ) goes like m%g\ X^/m%v^ A 2 <x A^/A 2 



6.2 Large Higgs portal coupling 

For larger values of A m , say A m > 10~ 3 , the h-rj mixing is larger and annihilations channels 
other than AA — > r\r\ in Fig. 1 become important, or even dominate the DM freezeout process. 
This leads to a more complex allowed parameter space. In agreement with the electroweak 
data and LEP constraints above, Fig. 4 displays sets of values of the parameters which lead 
to a relic density within the WMAP range above as well as to a direct detection cross section 
below the current upper limits [18] . We find such sets of values for h-rj mixing as large as 

4 This is similar to what happen's in "secluded" DM models [5], see also Ref. |26| . 
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Large A,„ regime 
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Figure 4: For 1CT 3 < A m < 1 and 114.4 GeV < m h < 180 GeV, values of g$ vs v^, rriA vs 
m v and vs cr(^4A — >• A/V) leading to 0.091 < O/i 2 < 0.129. A^ has been varied between 
10~ 5 and 1. Dots with ~ 2 mA proceed through resonance of the r\ exchange diagrams 
(to W + W~ , ZZ, ff or hh). Similarly dots with rrih ~ 2mA similarly are dominated by the 
Higgs exchange diagrams (to W + W~ , ZZ, ff or rjrj). Dots with m„ < 114.4 GeV are for 
suppressed values of sin /3 to agree with the LEP constraints on the h —> ff branching ratio. 



maximal. For such large mixing both h and r/ masses have to a very good approximation 
to be above 114.4 GeV to accommodate the LEP direct search limits and must remain low 
to accommodate electroweak precision data constraints, see section 5. In this case the LHC 
experiments should be able to discover both scalars, just in the same way as for the Higgs 
boson in the SM. As for the vector bosons they should also be presumably discoverable at 
LHC through virtual 77, although this will depend on their mass. Their mass can lie within a 
large range, even below the GeV scale in case the relic density can be obtained predominantly 
from the annihilation to fermions via h and/or ry exchange. 

As for the direct detection rate, we obtain a large range of cross section values, from 
several orders of magnitude below the current experimental limits to several orders of mag- 
nitude above them. Fig. 4 gives the sets of values of parameters which are consistent with 
these current limits. It would be worth to see if this model could also lead to interesting 
DM indirect detection signals (which, as for direct detection, are proportional to sin 2 /3 and 
therefore decouple from the relic density constraints for small mixing angle). In particular 
for what concerns positron indirect detection, as the 77 or h scalar mediated interaction be- 
tween 2 gauge bosons is attractive, it can lead to an enhancement of the positron emission 
through Sommerfeld effects [2TJ, [28] - It is to be analyzed in detail if, without applying large 
ad hoc boost factors, this enhancement can be large enough to explain the excess of positrons 
observed by the Pamela experiment [29]. From the results of Ref. |28| . Fig. 10 in particular, 
this appears to be possible, either for via in the multi TeV range through dominant anni- 
hilation to hh, hn, 7777, or for smaller values of tua (but larger than ~ 100 GeV) with m r] 
around ~ 1 GeV. In the latter case, rj particles (from A{Ai — > rjrj, hn annihilations) decay pre- 
dominantly to leptons (leading to a hard enough positron spectrum) and serve as mediators 
to produce highly non-relativistic leptons in large quantity, as in the mediator mechanism 
very recently proposed in Ref. [30j . Like in many models with such a light mediator this 
mechanism suffers nevertheless from stability problem of the light mediator mass under ra- 
diative corrections (i.e. in our case the stability of the 77 mass under gauge induced self-energy 



8 



corrections). 



7 Effects of higher dimensional operators 

The model above doesn't necessarily requires a UV completion^] A question one must ask 
nevertheless is, if there exist heavier particles in the hidden sector at a higher energy scale 
A, how these particles could spoil the stability of the vector bosons. If these particles induce 
dimension five operator destabilizing the vector bosons, one expects a lifetime, ta ~ cA 2 /m A 
(with c a coefficient of order 4tt for a two body decay) , many orders of magnitude smaller than 
the age of the universe, even for A as high as the GUT scale (unless the involved couplings 
are highly suppressed). However in the hidden vector model above there is no gauge invariant 
dimension 5 operators which could be induced by a higher energy physics. Only dimension 6 
operators can be induced, which induce much longer lifetimes, see e.g. [6]. 

For example one gauge invariant operator which could be induced is D ^ (j>D^ H / A 2 . 
Since (j> is assumed to be at a low scale this operator is relevant. It induces for example the 
decay A% — > (jxj)* with a lifetime of order ta ~ 47rA 4 /m^ which for m,A — 1 TeV (~ 1 GeV) is 
longer than the age of the universe if A is above ~ 10 13 GeV (10 9 GeV), or less if the involved 
couplings are smaller than unity. This means that the heavy particles at the origin of these 
operators must be heavier than these scales (just as for the proton in the SM in presence of any 
B violating new physics source, but at a much less constrained level because the lower limit 
on the proton lifetime is far larger than the age of the universe). Note that there are other 
dimension six operators which could cause the decay of the vector bosons, with similar lower 
bound on the underlying scale, ^F^f (f)F Y ^/A 2 , D^^D^H/A 2 , D u (f>F Y ^ /A 2 . 
There are also dimension 6 operators which do not violate the custodial symmetry, therefore 
not causing any decay, D^F'^ D u cf>/A 2 , D^^D^H^H/A 2 and F^F 1 ^ 'F Yu / A 2 Q 

It has been shown in Ref. [32] that DM decaying dominantly to e + e~ pairs with lifetime 
of order ~ 10 7 times the age of the universe, which we typically obtain from a dimension 
6 operator with A of order the GUT scale, leads to a positron flux with the right order of 
magnitude to explain the excess of positrons observed by the Pamela experiment. This leads 
also to a cosmic 7 flux that can be probed by current experiments Q 



5 In particular as it involves only one gauge interaction in the hidden sector and no fermions it doesn't call 
for any particular grand-unification UV completion in the hidden sector. As in the SM, there is though a 
hierarchy problem related to the stabilization of the mass of the Higgs boson 77 (under radiative corrections 
from any new physics or from the gravitation). We do not address this problem. 

6 Note also that if one adds in the hidden sector light fermions to which the hidden vector can decay, still 
one would get a DM candidate: the lightest fermionic multiplet of the non-abelian gauge-symmetry (and SM 
singlet) would provide a stable DM candidate because such hidden sector fermions can couple only in pairs to 
any other particle. This shows how generically gauged sectors coupling to the SM through the Higgs portal 
can provide DM candidates. It is beyond the scope of this work to study such possibilities which contain more 
fields than the hidden vector model, and are very different in many respects. 

7 Possible consequences of the decays induced by the operators above, in particular emission of intense 7 
lines, and excess of positrons, are analysed in Ref. |5T] . 



9 



8 A few more comments 



If one adds to the SM as few new fields as possible there are not that many possibilities to 
obtain a stable DM candidate without assuming by hand a discrete or global symmetry. 

One possibility, which holds with only one extra field [6], assumes the existence of a high 
fermion or scalar SU(2)l multiplet: a fermion quintuplet or higher or a scalar sextuplet or 
higher. Such multiplets are stable because no gauge invariant operators destabilizing these 
multiplets can be written with dimension less than 6. 

Another possibility, which holds with 2 extra fields, is the hidden vector model above. It 
involves lower multiplets. 

A third possibility, which holds with 3 extra fields, has been proposed in Ref. [9j. It 
assumes a U(l)' gauge boson, a scalar charged under the U{\)' and a fermion also charged 
under it. The scalar breaks the U{\)' to make the gauge boson massive and the fermion is 
the stable DM particle. The SM and hidden sector can communicate through both kinetic 
mixing and Higgs portal interactions]^] 

Finally note also that if in the hidden vector model above one considers an abelian gauge 
group G' = U(l)' (instead of the SU(2)hs above), with <fi a complex scalar charged under 
it (instead of the doublet above), one ends up with a Lagrangian similar to the one of 



Eq. (10), with nevertheless 3 important differences. First, by replacing the non-abelian field 



in Eq. (10) by the abelian one, the trilinear gauge couplings disappear and therefore all 
annihilation processes of Fig. 2 do not exist, but these processes are not mandatory to obtain 
the experimental relic density. Second with an abelian gauge symmetry there is no more 



custodial symmetry to make the gauge boson stable but still the Lagrangian of Eq. ( 10 ) with 
an abelian field instead of the non-abelian one displays a Z2 symmetry (under which is 
odd with all other fields even) related to the charge conjugation symmetry of the Lagrangian 
(which results from the gauge symmetry and particle content). Third, but not least, in order 
that the U{1)' gauge boson is stable one has to make the assumption that there is no F'FQ U 
kinetic mixing interaction with the hypercharge gauge boson. Unlike in the non-abelian case 
this term is not forbidden by any symmetry of the model. This cannot be justified without 
assuming extra symmetries. But it can be noted that, would this kinetic mixing term be 
absent, would all numerical results obtained above hold also for this case apart from factors 
of order unity (for instance a factor 1/3 in the relic density because there is only one DM 
component instead of 



9 Summary 

We have shown that a hidden sector vector multiplet associated to a non-abelian gauge group 
G', coupling to the SM only through the Higgs portal interaction of a scalar charged under 
this gauge group, constitutes a perfectly viable DM candidate. This vector multiplet is stable 
without needing to assume any discrete or global symmetry, due to the custodial symmetry 

8 Note that, with many more fields, mirror models which consider a complete copy of the SM in the hidden 
sector, can also lead to stable candidates (in particular the mirror proton) without assuming a discrete or 
global symmetry by hand, see e.g. Ref. [33] . 

9 The phenomenology of a similar model has been studied in Ref. [M] independently of DM considerations. 
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of the Lagrangian which results from the gauge symmetry. The stability of the gauge bosons 
could be spoiled by physics at higher energies but only through dimension six operators which 
is fine as long as the new physics is above a high scale, ~ 10 14 GeV (~ 10 10 GeV) for a DM 
mass equal to 1 TeV (1 GeV). For small Higgs portal coupling (but large enough to thermalize 
the hidden sector with the SM particles prior to freeze-out) the relic density is determined 
only by hidden sector parameters, while the direct detection rate necessarily involves the 
Higgs portal quartic coupling. This allows to decouple the direct detection constraints from 
the relic density ones. Large Higgs portal interaction, i.e. large mixing between the Higgs 
boson and the extra scalar, is also allowed in a large fraction of the parameter space, in 
case the model is testable at accelerators. The gauge boson mass can lie within a wide 
range of values from ~ MeV to few tens of TeV. This model can also lead to a Sommerfeld 
enhancement of the positron emission relevant for the recent Pamela experiment result. 
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A Appendix 

In terms of the various input parameters of Eq. Q, the parameters of Eq. ( 10 ) read as follows: 



m. 



K 



h 



K 



K h 

Si, 

& 

u 
\h 

M 
A 2 
A 3 

Py 
Ph 
Pi 
P2 



m\i cos 2 /3 + m 2 / sin 2 (3 — m h ,„, sin 2/3 
m 2 / cos 2 /3 + m\i sin 2 f3 + m 2 /^/ sin 2/3 

77^1 cos 2 (3 



; sin 2 /3 



-g0$ sin2 / 3 
^W/3 

^ 2 cos 2 /3 



1 



1 



-o 2 sin 2/3 

r 

;4 cos /3 

-g 2 sin (3 
g 2 cos /3 

A^, cos 4 /3 + A sin 4 /3 + A m cos 2 /3 sin 2 /3) 
A^, sin /3 + A cos 4 /3 + A m cos 2 /3 sin 2 /3) 

6A,* cos /3 sin 2 /3 + 6A sin 2 cos 2 /3 + A m (cos 4 /3 + sm 4 /3 - 4 cos 2 /3 sin 2 /3)) 

-4A,/, sin 3 /3 cos /3 + 4A cos 3 /3 sin /3 - A m sin 2/3 cos 2/3) 

-4A0 cos 3 /3 sin /3 + 4A sin 3 (3 cos /3 + A m sin 2/3 cos 2/3) 

4A0V0 cos 3 /3 + 4A-U sin 3 f3 + 2A m (t> cos 2 (3 sin /3 + sin 2 /3 cos /3)) 



-4A 



sin 



/3 + 4vA cos 3 /3 + 2A m (t> sin 2 /3 cos (3 — cos 2 /3 sin /3)) 



— 12u0A0COs /3sin/3 + 12Aw sin /3 cos (3 + A m (2t> cos /3 — 2v ( psm /3 — 4t> cos /3 sin /3 + 4^ cos /3sin/3) 
12A^U0 sin 2 (3 cos /3 + 12\v cos 2 /3 sin /3 + A m (2u sin 3 /3 + 2t>^ cos 3 /3 — 4v cos 2 /3 sin (3 — 4v,p sin 2 (3 cos /3)) 
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